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ABSTRACT 



Aims. Solar colors have been determined on the uvby-fS photometric system to test absolute solar fluxes, to examine colors predicted 
by model atmospheres as a function of stellar parameters (X e g, log g, [Fe/H]), and to probe zero-points of T eff and metallicity scales. 
Methods. New uvby-fi photometry is presented for 73 solar-twin candidates. Most stars of our sample have also been observed 
spectroscopically to obtain accurate stellar parameters. Using the stars that most closely resemble the Sun, and complementing our 
data with photometry available in the literature, the solar colors on the uvby-fi system have been inferred. Our solar colors are 
compared with synthetic solar colors computed from absolute solar spectra and from the latest Kurucz (ATLAS9) and MARCS model 
atmospheres. The zero-points of different T c ff and metallicity scales are verified and corrections are proposed. 

Results. Our solar colors are (b-y) e = 0.4105±0.0015, m liQ = 0.2122+0.0018, c UQ = 0.3319±0.0054, and/3 = 2.5915+0.0024. The 
(b-y)e and rai >0 colors obtained from absolute spectrophotometry of the Sun agree within 3-cr with the solar colors derived here when 
the photometric zero-points are determined from either the STIS HST observations of Vega or an ATLAS9 Vega model, but the Ci >Q and 
fl a synthetic colors inferred from absolute solar spectra agree with our solar colors only when the zero-points based on the ATLAS9 
model are adopted. The Kurucz solar model provides a better fit to our observations than the MARCS model. For photometric values 
computed from the Kurucz models, (b — y) Q and rai >0 are in excellent agreement with our solar colors independently of the adopted 
zero-points, but for c- lQ andyS agreement is found only when adopting the ATLAS9 zero-points. The c lG color computed from both 
the Kurucz and MARCS models is the most discrepant, probably revealing problems either with the models or observations in the u 
band. The r cff calibration of Alonso and collaborators has the poorest performance (~ 140 K off), while the relation of Casagrande and 
collaborators is the most accurate (within 10 K). We confirm that the Ramirez & Melendez uvby metallicity calibration, recommended 
by Arnadottir and collaborators to obtain [Fe/H] in F, G, and K dwarfs, needs a small (~10%) zero-point correction to place the stars 
and the Sun on the same metallicity scale. Finally, we confirm that the C\ index in solar analogs has a strong metallicity sensitivity. 

Key words, stars: atmospheres - stars: fundamental parameters - stars: solar-type - Sun: fundamental parameters 



1. Introduction 



Photo metry in the uvby-fi system dStromgrenl Il963t ICrawford 
1966) is well suited to the determination of basic stellar 
atmospheric parameters for F-, G-, and K-type stars through 
the color indices (b—y), m\ - (v-b) - (b-y) and c\ - (u- 
v) - {v-b). Several empirical calibrations exist in the lite rature 
to transform (b-y) o rffto ?W fe.g.lAlonso et aflll996L Il999t 
IClemetall l2004t Ra mirez & M elendez , l2005bi h ereafte r 
RM05b: iHolmberg et all l2007t ICasagrand eeiail l2010h . 
while the m \ index ca n be used to determine [Fe / H1 in 
dwarfs (e.g. [Stromgren, 1964; Gustafsson & Nissen, 1972; 



OlsenL \\9%4. ISchuster & NissenL 11 989; 



Malvuto, 1994; 
Haywood 1200 2; Martell & Laughlin, J200l Martell & Smith, 
2004 IRa mirez & Melendez, ~2005 al h ereafter RM05a; 
Twarogetall 120071: Holmberg et al„ 2007) 
Bond fl980t I Arellano Ferro & Mantegazzal 



3) and giants (e.g . 
1996; Hilker, 2000; 



* Based on observations collected at the H.L. Johnson 1.5 m tele- 
scope at the Observatorio Astronomico Nacional at San Pedro Martir, 
Baja California, Mexico. 



IRamirez & Melendez, 2004; Calamid a~et~aIL 120071 l2009h , as 
reviewed by lArnadottir et al.l (l2010t) . The evolutionary stage 
of stars can be determined us ing the ci index together with 
other uvby color indices (e.g. [Cr awford, 1975; Olsen, 1984; 
Niss en & Schusterl 11991b [Schuster etal., 2004; iTwarog et a~ . 
120071) . 

Since there are many difficulties in observing the Sun 
with the same in s trume n tation as w e obser v e other stars (e.g . 
Stebbins & KronL 1 19571 iGaUouetteL Il964t IClements & Neffl. 
1979; Tuegl 119821; lLockwood et all 1 1992b . the Sun cannot be 
used to set the zero-points of transformations between color in- 
dices and fundamental stellar parameters. Accurate transforma- 
tions are important in many areas of astrophysics. For exam- 
ple, in the study of the primordial lithium abundance, an accu- 
rate T e ff scale is needed to compare the Li abundance of metal- 
poor stars with that obtained from Big Bang Nuc leosynthesis 
(e.g. lMelendez et all l20loi ISbordone et all 1201 Oh . The terres- 
trial planet signatures found in the chemical composition of 
the S un (Melendez et al., 2009, her eafter M09: IRamirez et all 
2009, hereafter R09; Ramir ez et all l2010h shows that for accu- 
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rate comparisons between the Sun and stars an accurate temper- 
ature scale must be used in the determination of chemical abun- 
dances. Although there is inconclusive evidence about whether 
the Sun is too meta l-rich with respect to stars of simil ar age and 
Galactic orbit dHavwood 120081: iHolmberg et all l2009) implying 
that the Sun could ha ve been born in the inner part of the Galaxy 
dWielen et all 11996). this apparent offset between the Sun and 
stars could be partly due to zero-point errors in the photomet- 
ric metailicity scale (see e.g. Table 3 of Arnaddttir et al. 2010 
for systematic differences between spectroscopic and photomet- 
ric metallicities). Furthermore, we note that zero-point errors in 
the r e ff scale, as well as errors in the metailicity scale, would 
introduce systematic errors in the ages of stars determined from 
isochrones. 

According to their similarity to the Sun, stars can be 
classified as "solar-type stars" (late F to early K stars), "so- 
lar analogs" (G0-G5 dwarfs with solar metailicity within 
~ a factor of 2-3), and "so l ar tw i ns" (stars almost identi- 
cal to the Sun) (e . 
ISoderblom & Kingl 



Secchi, 1868; ICavrel de Strobe! 119961: 
19981: M09). Many works have used 
either solar-type stars or solar analogs to infer solar colors to 
improve or check the effective temperature scale, the perfor- 
mance of model atmospheres, and the absolute flux cal i bration 
of the Sun (e.g. [Petti t & Nicholso n. 1 19281: iKuiperl fl 938: 
Stebbins & Whitford, 1945; Stebbins & Kron, 1957; Johnson, 
1962; Kron, 1963; van den Bergh, 1965; Croft etal., 1972; 



OlsenL 1 1 976J ISchusterl 1 1 976t iBarrv et alTTl 97 8 : Har dorp1.il 978 



1980a b; Cle ments & Nefl . 
19811: iNe ckel & Labs! 
Tueg & Schmidt-Kaler, 1982 



1979 
1981 



o 



C96; 



2003 



Cavrel de Strobel et al 
Chmielewskil 1 1981 1 
198 31 iTavloil ll984J 
Vandenberg & Pol 



Mitchell & Schuster, 1985; Neckel, 1986; Vandenberg & Poll, 
1989 1: [GravL 119921 fl99iriStraizvs & Vahaugaf 1 19941 iTavlorl 
1994 ; ICavrel de Strobe iril996t IColina et all 1 19961 hereafter 



Cavrel de Strobel, 1996; Colina et al., 1996, hereafter 
Hauck & Kunzli, 1996; Bess ell et al.1119 98; Miro nov et al.[ 

20001: 
2006; 



1998; Sekiguchi & Fukugita 



RM05b; 



Melendez & Ramirez, 
Pasquini et all l2008t 



-.-_ . | H olmberg e t al ._ .. 
Rieke et all l2008t lCasagrand e et all 1201 Ob . 

Sola r twins h ave spectra almost identical to the Sun 
(Cavrel de Strobel, 1996), hence are better suited to setting the 
zero-points of fundamental calibrations, especially since they are 
found without assuming a priori a temperature scale, but their 
identification is based purely on a model-independent analysis 
of their spectra with respect to a solar spectrum obtained with 
the same instrumentation. However, until only a few years ago 
one solar twin was known (18 Sco; iPorto de Mello & da SilvaL 
1997; Soubiran & Tri audi 120041) . The situation has changed dra- 
matically with the ident ification of many st ars similar to the Sun 
(Melendez et al., 2006, hereafter M06; Melendez & Ramirez, 
2007, hereafter MR07; Takedaetal., 2007, hereafter T07; 



Takeda & Tajitsul 120091 hereafter T09; M09; R09; Baumann et 



al. 2010, hereafter B10). It is now feasible to use solar twins with 
accurately determined stellar parameters to test the predictions 
of model atmospheres and the accuracy of empirical photomet- 
ric calibrations. In the present work, we perform this study for 
the widely used uvby-B system. 



2. Photometric observations 

2. 1 . Selection of the sample 

Before starting our systematic survey of stars similar to the Sun 
in the Hipparcos catalogue, we performed pilot observations 
of solar-twin candidates to test our selection criteria. Our pilot 
study (M06) found zero-point offsets in our r e ff scale (RM05b). 



We therefore applied small corrections to the solar colors pre- 
dicted by RM05b to increase our chances of finding stars resem- 
bling the Sun. In particular, we searched for solar twins around 
the Tycho color (B - V) T = 0.7225 instead of (B - V) r =0.689 
predicted by our earlier color- T e s relations (RM05b). Our im- 
proved r e ff scale (Casagrande et al. 2010) indeed corrected the 
zero-point problem, and now predicts a solar (B - V);r=0.730, in 
good agreement with our tentative correction. 

The main parameters used to select solar- twin candidates 
from the Hipparcos catalogue were parallaxes and (B - V)t col- 
ors. Additional criteria used (when available) were other optical- 
infrared colors (e.g. Vj — K, b-y), photometric variability, infor- 
mation on multiplicity, previous literatu re values for [F e/H] (ob- 
tained from an updated version of the ICavrel de Strobel et all 
2001, catalogue), rotation, and chromospheric activity. We found 
initially about one hundred stars satisfying our selection criteria 
within 75 pes. We later expanded our search to cover the whole 
Hipparcos catalogue, increasing our sample by about 1/3. 



2.2. Observations and data reduction 

The uvby-B data presented here in Table 1 for the solar twins 
were taken using the H. L. Johnson 1.5 m telescope at the San 
Pedro Martir Observatory, Baja California, Mexico (hereafter 
SPM), and the same six-channel uvby-B photoelectric photome- 
ter as for the northern observations of Schuster & Nissen (1988, 
hereafter SN), for all the uvby-B observations of Schuster et al. 
(1993, hereafter SPC), the northern data of very-metal-poor stars 
by Schuster et al. (1996), the uvby-B data for very-metal-poor 
stars in Table 1 of Schuster et al. (2004), and the uvby-B data 
for high-velocity and metal-poor stars in Table 1 of Schuster et 
al. (2006). The new uvby-B values for solar twins included here 
in Table 1 were taken during three observing runs in November 
2007 (8 nights), April 2008 (7 nights), and September 2008 (4 
nights). 

The uvby-B solar-twin data were taken and reduced using 
techniques very nearly the same as for SN and SPC (see these 
previous papers for more details). The four-channel uvby sec- 
tion of the SPM photometer is really a spectrograph-photometer 
that employs exit slots and optical interference filters to de- 
fine the bandpasses. The grating angle of this spectrograph- 
photometer was calibrated using a cadmium lamp at the begin- 
ning of each observing run to position the spectra on the exit 
slots to within about +1A. Whenever possible, extinction-star 
observations were made nightly over an air-mass range of at least 
0.8 (see Schuster & Parrao 2001; also Schuster et al. 2002), and 
spaced throughout each night several "drift" stars were observed 
symmetrically with respect to the local meridian (two hours east 
and then two hours west). Using these observations, the atmo- 
spheric extinction coefficients and time dependences of the night 
corrections could be obtained for each of the nights of observa- 
tion (see Gr0nbech et al. 1976). Finding charts were employed 
at SPM to confirm identifications of the program stars and to 
select regions for the "sky" measurements. As for previous stud- 
ies, such as SN and SPC, the program stars were observed at 
SPM to at least 50 000 counts in all four channels of uvby and 
to at least 30 000 counts for the two channels centered at HB. 
For all program stars, the sky background was measured until 
its contributing error was equal to or smaller than the error in 
the stellar count. At SPM, an attempt was made to obtain three 
or more independent uvby observations for each of the program 
stars, i.e. photometric observations during at least three inde- 
pendent nights; this aim was achieved, or exceeded, for all solar 
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twins except HIP75923 and HIP77883 for which we obtained 
two observations each. 

As for the SN and SPC catalogues, all data reduction was 
carried out following the precepts of Gr0nbech et al. (1976) 
using computer programs kindly loaned by T. Andersen (see 
Parrao et al. 1988). At SPM, the uvby-fi standard stars observed 
were taken from the same lists as for the previous catalogues, 
and are mostly secondary standards from the catalogues of Olsen 
(1983, 1984). The reduction programs create a single instru- 
mental photometric system for each observing run, including 
nightly atmospheric extinctions and night corrections with lin- 
ear time dependences. Then transformation equations from the 
instrumental to the standard systems of V, (b-y), m\, c\, and/? 
are obtained using all standard stars observed during that ob- 
serving period. The equations for the transformation to the stan- 
dard uvby-fi system are the linear ones of Crawford & Barnes 
(1970) and Crawford & Mander (1966). Small linear terms in 
(b-y) are included in the standard transformation equations for 
m\ and c\ to correct for bandwidth effects in the v filter. Our y 
measurements were transformed onto the V system of Johnson 
etal. (1966). 

Thirty-six uvby-fi standard stars were employed during the 
observing run of November 2007 providing instrumental photo- 
metric errors ranging from 0.002 mag. in (b-y) to 0.009 mag. in 
c\, and errors in the transformations to the standard photometric 
system from 0.004 to 0.01 1 magnitudes, respectively. For April 
2008, these values were for 33 standard stars, 0.002-0.010 mag., 
and 0.005-0.013 mag., respectively, and for September 2008, for 
35 standard stars, 0.002-0.008 mag., and 0.006-0.009 mag., re- 
spectively. Instrumental and transformation errors in the magni- 
tude V and in the indices m\ and/3 were always intermediate in 
value between those given above. 



2.3. The catalogue 

Table 1 presents the uvby-fi catalogue for the 73 solar-twin can- 
didates observed at SPM. Column 1 lists the Hipparcos number; 
Col. 2 gives the HD (or BD) number, Col. 3 the V magnitude 
on the standard Johnson UBV system; and Cols. 4-6 and 9, the 
indices (b-y), mi, c\ and fi on the standard systems of Olsen 
(1983, 1984), which are essentially the systems of Crawford & 
Barnes (1970) and Crawford & Mander (1966) with north-south 
systematic differences corrected. Columns 7, 8, and 10 give Ny, 
N uv b y , and Np, the total numbers of independent V, uvby, and fi 
observations. 

A very small subset of our photometric observations was 
made through light cirrus clouds in the absence of moonlight. 
It has been well documented (e.g. SN; Olsen, 1983) that obser- 
vations in the indices b-y, ni\,c\, and fi made with simultaneous 
multichannel photometers are not affected in any significant way 
by light (or even moderate) cirrus, while the V magnitude, ob- 
tained from only the y band, is affected. For this reason, a few of 
the solar twins, such as HIP60314, HIP74389, and HIP1 18159, 
have fewer independent observations of the V magnitude than 
the indices. 



2.4. Comparisons with other photometric data 

We searched for uvby-fi photometry in the G eneral 
Catalogue of Photometric Data (iMermilli od et all Il997t 
lHauck & Mermilliodl 119981) and found data for 23 solar-twin 
candidates and 12 solar analogs (observed by ourselves spectro- 



scopically for other projects). The photometry of these 35 stars 
is given in Table 2. 

The accuracy of uvby-fi photo metry obtained at SPM has 
been extensiv e ly tes t ed (e. g. SN; lArellano Ferro et all [1990; 
ISchuster et~ail I1993L 119961) . and we illustrate below that our 
solar-twin photometry is also in excellent agreement with the lit- 
erature, in all cases with mean differences well below 0.01 mag- 
nitudes. There are 12 stars in common between our sample and 
previous work. The average difference (ours - literature) in (b-y) 
is only -0.001 (<x = 0.004). Our V Johnson photometry is also 
in good agreement, with a mean difference (ours - literature) of 
only +0.001 (cr = 0.006). The colors m\ and c\ also compare 
well, with a difference (ours - literature) of only +0.002 (cr - 
0.006) and +0.006 (cr = 0.015). 

Considering the excellent agreement between the photom- 
etry available in the literature and in our own data sets, and 
considering the previously shown accuracy and precision of the 
SPM uvby-fi photometry, we conclude that both data sets, of 
Tables 1 and 2, provide solar-color indices very close to the stan- 
dard uvby-fi system. 



3. The solar-color indices 

Our observations comprise the largest photometric data set yet 
taken of stars very similar to the Sun in the uvby-fi system. 
Adding other photometry available in the liter ature and using 
our ow n accurate stellar parameters (M09; R09; Melendez et all 
2010b, hereafter Ml 0b), and taking into account the variations 
in colors with T e g, log g and [Fe/H], we can infer the "solar" 
colors by interpolating them to the stellar parameters of the Sun : 
7W = 5777 K, log g = 4.44, and [Fe/H] = 0.00 (e.g. lCoxl 120001 
Gray, 2005; RM05b). 

The quality of our stellar parameters (M06; MR07; M09; 
R09; MlOb) is very high because both the Sun (reflected light 
of asteroids) and the solar twins were observed with the same 
instrumentation during the same observing runs, and all data 
reduction and analysis were performed in an identical way. 
Our spectra have typically a resolution of 60,000 and S/N~200 
for stars observed with the 2.7m telescope at McDonald and 
S/N~450 for stars observed with 6.5m Magellan Clay telescope 
at Las Campanas. Errors as low as ~25 K in r e ff, 0.04 dex in 
log g, and 0.025 dex in [Fe/H] can be obtained in the best cases, 
and abundance ratios with errors as low as 0.01-0.02 dex have 
been obtained with the above data, showing that the Sun is a star 
with a peculiar chemical composition (M09; R09). Additional 
spectra available in the literature for solar-twin candidates were 
analyzed by B10 using similar techniques. The spectra are from 
HARPS obser vations available at th e ESO archive and from the 
S 4 N database dAllende Prieto et al .. 2004fl In both cases, a so- 
lar spectrum taken with the same instrumentation was employed 
in the differential analysis. 

The stellar parameters were determined homogeneously by 
our own group (M06, MR07, M09, R09, Ml 0b, B10) using 
Kurucz model atmospheres and a line-by-line differential anal- 
ysis with respect to the Sun. The r e ff, log g, [Fe/H], and micro- 
turbulence were determined iteratively until both the differential 
excitation equilibrium of Fel lines and the differential ioniza- 
tion balance of Fel and Fell, were achieved. The microturbu- 
lence was also determined simultaneously, by requiring no de- 
pendence of the iron abundance (from Fel lines) on the reduced 
equivalent width. 



1 S 4 N: Spectroscopic Survey of Stars in the Solar Neighborhood; 
available online at http://hebe.as.utexas.edu/s4n/ 
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Since our solar-twin sample spans a relatively narrow range 
in atmospheric parameters relative to the Sun, in principle sim- 
ple linear fits of color versus (vs.) each parameter (T e ff, log 
g, and [Fe/H]) would provide sufficiently good estimates of 
the solar colors. As suggested by the referee, a global fit (e.g. 
Mitchell & Schuster, 1985) to all stellar parameters (color = 
f(T e s, log g, [Fe/H])) would be preferable due to the mutual in- 
terdependence of the stellar parameters. Fortunately our sample 
includes also solar analogs covering a broader range in colors 
and stellar parameters than the solar twins, so that the depen- 
dence on the different stellar parameters can be well determined 
by a global fit. The following formula was employed: 

color = A + B (T eS - 5777) + C (logg - 4.44) + D [Fe/H]. (1) 

The advantage of this equation is that A will give us directly the 
solar color, while its uncertainty could be determined from the 
error in A or from the scatter of the fit. 

From our sample of stars, we selected a group of solar twins 
with r eff , log g, and [Fe/H] within 100 K, 0.1 dex, and 0.1 dex 
of the solar parameter^ given above. To perform a more robust 
global fit, we extended our solar twin sample with solar analogs 
covering a range in T e g, log g, and [Fe/H]. Thus, to explore the 
metallicity dependence of the colors, we selected a sample of 
solar analogs with the same constraints in T e ff and log g, but 
covering a broader range in metallicity (-0.4 < [Fe/H] < -0.1 
dex, +0.1 < [Fe/H] < +0.4 dex). In a similar way, to improve 
the fit to T e g, we selected a group of solar analogs with the same 
constraints on log g and [Fe/H] as the solar twins, but with T e ff 
in the range (5617 < T eS < 5677 K, 5877 < T eff < 5937 K). 
Finally, to fit the trend with log g, we used a sample of solar 
analogs with the same constraints on T e ff and [Fe/H] as the solar 
twins, but covering a broader range in log g (4.29 < log g < 4.34 
dex, 4.54 < log g < 4.59 dex). 

The stellar parameters were taken from our work on solar 
twins (M06; MR07; M09; R09; MlOb; B10) and complemented 
in so me cases with other accu rate values availa ble in the litera- 
ture jlValenti & FischerL I2005L hereafter VF05; ILuck & HeiterL 
2006, hereafter LH06; T07: ISousa elaH 120081 hereafter S08; 
T09). The adopted stellar parameters for the solar twins and so- 
lar analogs are given in Table 3. 

The global fits to (b-y) had only three outliers, HIP 7245, 
HIP 8 15 12, and HIP 88427 which seem too red in (b-y ) for their 
T e ff and [Fe/H]. Using the lKaratas & Schuster! (12010) intrinsic- 
color calibration, we find that these three stars may be slightly 
reddened (E(b-y) ~ 0.015 - 0.023), although we have to bear 
in mind that the accuracy of the reddening calibration is of the 
same order. These stars were removed from the global fits. The 
results from the fits are presented in Table 4 and in Figs. 1-4, 
where filled circles represent the solar twins and open circles the 
solar analogs across a broader range of stellar parameters. 

3.1. The (b-y) solar color 

The global fit of (b-y) (Fig. 1 and Table 4) shows strong depen- 
dences on r e ff (at the 14-cr level) and [Fe/H] (10-cr). There is no 



2 Strictly speaking, a solar-twin star must be identical to the Sun 
within the observational uncertainties. Our definition based on derived 
stellar parameters is simply more practical because it is less dependent 
on the quality of the stellar and solar spectra available. 

3 The star HIP 79186 also seems reddened according to the 
Karatas & Schuster ( 20 1 ) calibration, but the stellar parameters of this 
star (r efl , log g, and [Fe/H]) = (5709 K, 4.27 dex, -0.12 dex) (R09) do 
not fall either in our solar-twin or our solar-analog samples, so it was 
not considered in the global fits. 
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Fig. 1. (b-y) vs. T e fi (upper left panel), [Fe/H] (upper right 
panel), and V magnitude (lower left panel, with a dot-dashed 
line at (b - y) Q ). The results of the global fit vs. the observed 
(b-y) color is presented in the lower right panel, with the dotted 
line indicating equality. Solar twins and solar analogs are repre- 
sented by filled and open circles, respectively. The dependences 
of the fit on T^ and [Fe/H] are shown by dashed lines, while 
the relative predictions of MARCS models (normalized to our 
inferred solar colors) are shown by solid lines. 



dependence (within the errors) on log g, therefore this parameter 
was excluded from the global fit. The star-to-star scatter from 
the fit is only 0.005 magnitudes, which is what is expected from 
the observational uncertainties (0.004-0.006 magnitudes). 

The standard error (s.e.) in the solar color was obtained from 
the observed star-to-star scatter and the number of data-points 
(s.e. = cr/ -^/sample size). The error in (b-y) B (and the other solar 
colors) was conservatively adopted as three times the standard 
error. Thus, we propose 



(b-y) = 0.4105 (±0.0015). 



(2) 



A plot of (b-y) vs. V magnitude can help us to reveal whether 
fainter (i.e., more distant) stars bias the above derived solar color, 
due to possible interstellar reddening. This plot is shown in the 
lower left panel of Fig. 1 . As can be seen, there is no trend with V 
magnitude. A linear fit of (b-y) vs. V, indeed shows a zero slope 
within the errors (slope = 0.0003 ± 0.0008). We note that most 
stars of the sample used in the global fit are brighter than V = 9, 
i.e. closer than ~ 68 pc. Even the faintest stars (V ~ 9.3) extend 
only to ~78 pc. According to Nal interstelar absorption maps, 
very little Nal abso rption is detected for distances up to ~ 80pc 
from the Sun dLallement et all l200l IWelsh et all l20loh . Thus, 
most of our sample is not expected to be significantly affected 
by interstellar absorption. As already mentioned above, the few 
stars that show some small sign of interstellar reddening were 
not included in the global fits. 
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Fig. 2. j3 vs. r e jf (upper left panel), [Fe/H] (upper right panel), 
and V magnitude (lower left panel). The results of the global fit 
vs. the observed /3 color is presented in the lower right panel. 
Symbols are as in Fig. 1, except that two curves are presented 
for MARCS models: the shallower curve is the one computed 
by us while the other one is from O09. Neither curve is able to 
explain the low/? colors observed for [Fe/H] < -0.2. 

3.2. The (3 solar color 

The global fit of/3 shows a dependence on r e ff, log g, and [Fe/H] 
(Table 4). Interestingly, the strongest dependence is with [Fe/H] 
(slope significant at the 5-<x level), while the dependence with 
both r e ff (2.9 cr) and log g (2.7 cr) is only at the 3-<x level. In 
Fig. 2, we show the dependence of the ft color index with respect 
to r e ff and [Fe/H]. The scatter from the fit is 0.006 mags, which 
is compatible with the observational errors in j3. 
The adopted ft solar color is 



Pq = 2.5915 (+0.0024). 



(3) 



Although there is a clear trend of /3 vs. V magnitude (Fig. 
2), it is not related to reddening because the [5 index is not af- 
fected by interstellar absorption. The bright stars (V < 6) that 
are causing the trend are within ~17 pc, hence they are not red- 
dened. The bright (V < 6) solar twins and analogs falling sys- 
tematically above the derived y6 , are stars hotter than the Sun, 
or more metal-rich than the Sun (or both), and therefore with ayS 
color systematically higher than solar because [5 increases with 
both increasing metallicity and T e ff. We note that the trend seen 
in nearby stars for f3 is not present in (b-y) because this color 
has opposite trends with T e ^ and [Fe/H]. Thus, by performing a 
deeper solar-twin survey than previous works, we avoided sys- 
tematic biases (such as selecting mainly hotter and more metal- 
rich stars) that might have been present if only brighter stars 
were studied. 



3.3. The mi solar color 

It is well known that the ;«i color correlates very well with 
metallicity (see references in the introduction), and this is clearly 
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and V magnitude (lower left panel). The results of the global fit 
vs. the observed m\ color is presented in the lower right panel. 
Symbols are as in Fig. 1. 



shown in Table 4, where according to the global fit the depen- 
dence on [Fe/H] is significant at the 21-cr level. The second most 
important parameter is T^ ( 1 cr), but log g also produces an im- 
portant dependence (slope significant at the 5-cr level). 

The tight correlation between m\ and [Fe/H] is shown in 
Fig. 3. The star-to-star scatter from the global fit is only 0.006 
magnitudes, which is compatible with the observational uncer- 
tainties. The m\ solar color recommended for the Sun is 

wi, = 0.2122 (+0.0018). (4) 



3.4. The c\ solar color 

In Fig. 4, we show the relation between c\ and log g, which is 
signi ficant at the 4 - cr level. As already n oticed in the literature 
(e.g. iTwarog et all l2002t IQnehag et all I2009L hereafter O09), 
the c\ index in late G dwarfs has a sensitivity to metallicity. Our 
global fit (Table 4) confirms the dependence on [Fe/H], which 
is actually more significant (8-<x) than the dependence on log g 
(Fig. 4). 

The predicted c\ solar color from the global fit is 



ci, Q = 0.3319 (+0.0054). 



(5) 



The star-to-star scatter in the global fit is 0.0165, which is con- 
siderably larger than the observational errors (0.009 - 0.013) 
for the c\ index. Following the suggestions of the referee, we 
explored whether this index is particularly sensitive to either 
anomalies in the chemical composition or the microturbulence 
velocities. 

Previous works show the effects of C and N upon the 
Stromgren 4-color uvby filter-measurements, via the NH and CN 
bands. For example, it seems that the NH band at 3360A affects 
the u measurements while CN affects the v measurements. Thus, 
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the c\ color that depends o n both u and v should be a f fected 
by abundance anomalies (Bond&Neff, 1969; Bond, 1974; 
Zacs et al l 1998; Grundahlet al., 2000, 2002; Schuster et al., 



2004: lYong et all l2008t) . In this context, it would be important 



to assess whether the small abundance anomalies in the solar 
chemical composition (M09; R09), in particular the difference 
between the highly volatile elements ( C, N, O) and Fe, may af - 
fect the uvby colors. As discussed by Strom gr en et alj 0982), 
variations in the He abundance may also affect the c\ index. 
Although the Hyades c\ anomaly for stars with (b - y) close 
to solar was found initially to be Acj ~0. 03-0. 04 (Hyades - 
field stars, or Hyades - Coma), it seems that the anomaly may 
only amo unt to Aci =0.024-0.02 5 after instrumental effects are 
corrected dJoner & Tavloii Il995l) . Another i mportant parameter 
affecting the colors may be microturbulence dConti & Deutsch, 
1966; NissenL ll98ll) . 

To test the above effects, synthetic spectra were computed for 
solar twins with variations in A[C,N,O/Fe]=-0.05 dex (M09), 
Av m j cro = +0.1 km s (most solar twins and close solar analogs 
have Vmicro within +0.1 km s _1 of the solar value), and an in- 
crease of 10% in the He abundance ( by number). Fluxes were 
computed with the code SYNTHE dKurucz & Avrettl Il98ll 
ISbordone et"ail 12004 iKu rucz, 2005) usin g ATLA S 12 model at- 
mospheres dKuruczL 1 19961 120051: ICastelht 120051) computed for 
the different aforementioned assumptions. The atomic line list 
adopted i n the spectr al computati ons is based on the com pila- 
tions bv lCoelhoet"ai] d2005l) and lCastelli & Hubrigl d2004 . and 
the mole cules C ? CH, C N, CO, H 2 , MgH, NH, OH, SiH, and 
SiO from Kurucz ( 1993) were included. The change in He does 
not significantly affect the c\ index (<0.001 magnitudes), but the 
change in microturbulence increases c\ by +0.0035 mag., while 
the change in C,N,0 increases c\ by +0.003 magQ Thus, small 
changes in chemical composition and microturbulence may help 
to explain the extra scatter seen in the global fit. Casagrande et 
al. (in preparation) demonstrate that the Stromgren system is not 
only sensitive to [Fe/H], but that it is possible to obtain infor- 
mation about the [a/Fe] ratio using uvby photometry. The Aci 
anomaly in the Hyades could be due to the effect of metallic- 
ity on c\. According to Friel & Boesgaard (1992), Hyades has 
an iron abundance 0.18 dex higher than Coma, which accord- 
ing to Table 4 corresponds to Acj = +0.022 mag., which is 
very cl ose to the Hyades anom aly relative to Coma (0.024-0.025 
mag.) dJoner&Tavlorlll995l) . 



4. Comparison with previous empirical 
determinations 

As can be seen in Table 5 (lower part), our inferred (b-y), m\, c\, 
and/3 solar colors agree well (within l-<x) with most earlier em- 
pirical results. However, previous determinations of solar col- 
ors have much lower control and homogeneity in T e g, log g, 
and [Fe/H], thus in some cases the agreement (within the error 
bars) could be due to fortuitous compensations of T e g, log g, and 
[Fe/H] effects on colors, or due to large errors in other studies. 

Several previous works (Saxner & Hammarback, 1985; 
iGravi 119921: Edvardsson et all 119931 Casagrande et al., 2010) 
give (b—y)o values that agree with ours within 0.005 m agnitudes, 
but o t her determinations are either too blue dCav rel de Strobel, 
[19961 IClem et all 2004; RM05b; IHolmberg et all 120061) or 
too red (Gehren, 1981). A previous empirical result for /? 
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Fig. 4. c\ vs. log g (upper left panel), [Fe/H] (upper right panel) 
and V magnitude (lower left panel). The results of the global fit 
vs. the observed c\ color is presented in the lower right panel. 
Symbols are as in Fig. 1. 



dSaxner & Hammarbacld.ri985L n = 2.591 ± 0.005) is in excel- 
lent agreement with our valu e (2.5915+0.0024). Previous deter- 
minations of wi >0 and ci Q (Cle m et all 120041: IHolmberg et all 
2006) also agr ee with o urs within the errors, but note that the 
errors given by IHolmberg et al] d2006l) are very large (0.03 and 
0.07 for ;«i and c\, respectively), so their solar colors cannot 
be used to perform stringent tests of model atmospheres and the 
temperature and metallicity scales. 

Although our inferred solar colors agree with previous de- 
terminations based on solar-type stars, our values based on so- 
lar twins should be preferred because of their more accurate 
and precise stellar parameters, resulting in correspondingly ac- 
curate and precise inferred solar colors. We use our "solar" col- 
ors below to test absolutely calibrated solar spectra, model at- 
mospheres and the r e ff and metallicity uvby-fi scales. 



5. Comparison with measurements on reflected 
(asteroid) solar spectra 



4 The changes in (b - y) and n\\ due to changes in He, CNO, and 
microturbulence are even smaller than for C\ . 



lOlsenl (1 19761) measured the J3 Q color of reflected solar light from 
asteroids, which seems justifiable because the /3 index should be 
largely independent of the wavelength dependence of the albedo. 
His measurement (J3 Q - 2.5955 ± 0.0024) is in good agreement 
(withi n 1-cr) with our inferred value (J3 = 2.5915 ± 0.0024). 
lOlsenl dl976l) did not measure the solar (b-y), but indirectly 
estimated (b - y) G from a color transformation using f3 & . His 
(b - y) Q = 0.390 ± 0.004 disagrees with ours (0.4105+0.0015), 
but this is probably due to the errors introduced by his adopted 
transformation fro m (3 to (b-y). Using the (b - y)o~/3 intrinsic 
color calibration of lKaratas & Schuster! d2010l) . the/3 = 2.5955 
from Olsen ( 1976) indeed implies that (b-y) e = 0.4062, in good 
agreement with our (b-y) . Using the same transformation, our 
/Jo = 2.5915 gives (b-y) & = 0.4089. 
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Unfortunately, direct measurements of the (b-y), m\ and c\ 
color indices for reflected solar spectra are not very useful be- 
cause of the color of the asteroid albedosQ A comparison with 
solar colors inferred from absolute spectrophotometry of the Sun 
is presented below. 



6. Comparison with synthetic colors 

Our accurate and precise colors inferred for the Sun can be used 
to test the performance of theoretical solar flux models and the 
quality of absolute solar flux measurements. However, additional 
ingredients enter the computation of synthetic colors, namely the 
adopted set of filters and the flux/magnitude zero-points adopted. 
We discuss them in the following. 

Since intermediate-band Stromgren filters are centered on 
specific spectral features, a correct characterization of the to- 
tal throug hout becomes crucial for generating synthetic colors 
(e.g. Lester et all Il986t O09). To test the influence of the fil- 
ter transmission curves on our results, we computed uvby in- 
dices using two different se ts of pass-bands, the original ones of 
ICrawford & Barnes! d 19701) and a s et that should be more repre- 
sentative of the SPM observations dBessellll2005l) . 

The p index is defined as the ratio of the flux measured 
through narrow (half- width of about 30 A) and wide (about 
150 A) profiles both centered on the H/3 line. In this case, 
the (212,214) filter transmission curves (Crawford & Mander, 
1966), the photomultiplier sensitivity, the atmos pheric trans mis- 
sion, a nd the reflectivity of aluminum given in Castelli & Kurucz 
(2006D were used to generate the total Tp throughout, ac- 
cording to the prescriptions of the beta.forcd program at the 
Kurucz websitqj. Indices calculated using Tp define the nat- 
ural system /?', which should be transformed using a set of 
equations to ag ree to the standard s ystem yS defined by the 
observations of ICrawford & Mander (1966). For the filter set 
(212,214), the transformatio n equation isyS = 0.248 + 1.368/3' 
( ICrawford & Manderl[l966h . 

In principle, to mimic the uvby-/3 observations presented 
here, synthetic photometry should reproduce the SPM instru- 
mental system, and the same transformation equations should 
then be applied to generate the standard system. In practice, 
this can hardly be done, since the SPM instrument is a six- 
channel spectrophotometer, and we follow the approach nor- 
mally adopted in the literature, i.e. of reproducing the standard 
system directly, by fixing the zeropoints using Vega. In Sect. 
2.4, we have shown the excellent agreement between our ob- 
servations and other photometric measurements, meaning that 
the transformation from the instrumental to the standard system 
is indeed accurate; we therefore expect our approach to return 
meaningful synthetic colors. 

The spectral energy distribution adopted for Vega and its ob- 
served indices also affect the outcome of synthetic photometry 
in the proce ss of establishing the ze ropoints of a photometric 
system (e.g. Casagrande et al., 2006). We used a spectrum ob- 
tained with the STIS spectrograph onboard the HST ([Bohlin, 
2007) of resolution 'R - 500, which enables the highest accu- 
racy (~ 1%) measurements achievable to date, and adopted the 
following averaged values for Veg a: (b-y) — 0.003, m\ —0.157 , 
d = 1.088 and/3 = 2.904 from lHauck& Mermilliodl (Il998l) . 



5 Although the colors of the Sun measured using asteroids may be 
affected, our spectroscopic analysis based on high-resolution spectra 
should not because we measure the flux relative to the adjacent contin- 
uum for narrow spectral lines 

6 http://kurucz.harvard.edu/programs/COLORS 



If we had adopted the colors of Vega given in ICrawford et alJ 
(Il972l) . the differences would have amounted to 0.001 magni- 
tude at most. 

Despite the complications posed by the pole-on and rapidly 
rotating nature of Vega, the effects on the blue part of 
the spectrum are expected to be sm all or negligible (e.g. 
ICasagrande et all [2006; Bohlin, 2007, and references therein). 
It is more relevant that we use the observed HST spectropho- 
tometry of Vega; additional comments on this issue are made in 
the following subsections. 

6.1. Colors from Kurucz and MARCS models 

Although synthetic uvby-fi colors computed using earlier 
MARCS and Kurucz models are found in the literature (see 
Table 5), we feel it appropriate to make our comparisons us- 
ing the most recent rele ases available (Castelli & Kurucz, 2004; 
iGustafsson et all I2008I) . In addition, as we have already men- 
tioned, different ingredients enter the computation of synthetic 
colors, and we differ from most of the previous works in that we 
use an observed spectrum of Vega to define the zeropoints. In 
principle, this should be the best approach for replicating obser- 
vations and for a correct comparison of solar flux models and 
absolute measurements with the colors determined from our so- 
lar twins. This choice mimics the observational approach, and 
the successes or failures of synthetic colors depend mostly on 
the quality of the solar input spectra. In practice, the situation is 
less clear, as we discuss further below and in Sect. I6.2I 

When computing synthetic colors from theoretical models, 
the use of a model atmosphere to describe also Vega may have 
the advantage of (partly) compensating for model inaccuracies 
by including these in the zero-points. 

We performed this exercise by taking from the Kurucz web- 
site the latest ATLAS9 model fluxes for Vega and the Sun; be- 
cause of the internal consistency of this approach, it should also 
determine which set of filters should be used. For the colors of 
the Sun we obta in (b-y) - 0.413, nu = 0.2 36 and c\ = 0.297 us- 
ing the filters of lCrawford & Barnesl (fl9 70). and (b-y) = 0.406, 
mi = 0.214, c\ = 0.303 for the Bessell (200J) filters. While (b- 
v) is reproduced with the two set of filters, the latter set of fil- 
ters provides results that are more comparable to our measured 
m\ B and ci , and therefore in the following we consider only 
the Bessell (2005) passbands. The effect of instead using those 
of lCrawford & Barnesl (1 19701) can be easily estimated from the 
above differences. 

The resolution of the spectra used to generate synthetic col- 
ors may also in principle affect the results. We tested ATLAS9 
model fluxes for Vega and the Sun at various resolutions ranging 
from % - 500000 to 200 and verified that for various combina- 
tions of these, differences in (b-y), mi and c\ always lie below 
0.001 magnitude]] 

However, for "R < 2000 the synthetic /3 index scales differ- 
ently for Vega and the Sun, i.e. even if synthetic spectra of the 
same resolution are used for the two stars, the value obtained 
for ft depends on 'R to an extent that may vary from a few milli- 
magnitudes up to several hundredths of a magnitude. For exam- 
ple, using high resolution ATLAS9 spectra for both Vega and the 
Sun, we obtain j6 = 2.587, which is in good agreement with our 
solar value. Using instead synthetic spectra of Vega and the Sun 



7 The Castelli & Kurucz (2004) grid of fluxes has a resolution vary- 
ing from 150 to 250 in the wavelength region of interest, and for this rea- 
son we took instead model fluxes at higher resolution from the Kurucz 
website 
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at % — 500 sampled at the same wavelength points, we obtain 
ft - 2.601 from which we estimate a (model-dependent) correc- 
tion of 0.014 magnitude when working at this low K. 

Because of this limitation on the STIS resolution of the 
Balmer line, the /3 indices computed for Table 5 have been ob- 
tained as follows: the solar spectra were downgraded to % — 500, 
sampled at the same wavelength points as the STIS spectrum, 
and the aforementioned correction was then applied. Using in- 
stead a high resolution ATLAS9 spectrum of Vega to define the 
zero-points does not require a downgrade to the resolution of 
the solar spectra, and the /? indices are changed by -0.030 mag- 
nitude (i.e. are bluer). We verified that for the other Stromgren 
indices the STIS resolution is high enough. The effect of using 
an ATLAS9 model flux of Vega to set the zero-points amounts 
to -0.004, 0.006, and 0.027 magnitudes for the (b-y), m u and c\ 
indices, respectively. In Figure [5] synthetic colors obtained with 
these two choices for Vega, i.e. STIS vs. ATLAS9, are compared. 

If the STIS spectrum of Vega is used, the Kurucz solar flux 
returns a value of (b-y) in excellent agreement with our ob- 
served value, whereas a considerably bluer color is obtained with 
MARCS. The greater accuracy of the Kurucz with respect to the 
MARCS solar flux in b and y bands can be n oticed also from th e 
spectrophotometric comparison presented in Edvardsson (2008). 
For mi, the Kurucz model is still in better agreement than the 
MARCS, one being slightly bluer and the other redder with re- 
spect to our mi . The same bluer and redder performance is 
also obtained for c\, although both models provide a rather poor 
match. The /3 index is similar for both models and considerably 
redder than our/? . 

In general, the ATLAS9 solar flux model performs better 
than the MARCS model, but they both have problems in repro- 
ducing ci jQ (in opposite directions) andyS (both systematically 
redder). Changing to the ATLAS9 flux of Vega to set the zero- 
points has a negligible impact on (b-y) and mi, but brings the 
theoretical /3 index in almost perfect agreement in both cases. 
For the c\ index, only the ATLAS9 result is helped by this 
choice, and this could partly stem from compensating errors in 
the ATLAS9 models of Vega and the Sun. 

In Table 5, we also show the color indices computed by O09 
using MA RCS models (also for the flux of Vega) and the pass- 
bands of lCrawford & Barnes! (Il970h . Their (b-y) agrees with our 
synthetic one, but their m\ and c\ colors are somewhat redder 
and bluer, respectively, in a way that is consistent with the dif- 
ferent passbands they adopted. The c\ solar color computed by 
O09 is in excellent agreement (within 1-cr) with our ci . 

6.2. Colors from solar spectra 

We computed (b-y), ni\, c\, and /3 indices using the absolute 
measurem ents of solar sp e ctra b y C96, Neckel (199 9, hereafter 
N99) and iThuillier et al.1 d2004l) . The spectrum of lRieke etaT] 
(2008) was not employed because of its very low resolution (R ~ 
100). The STIS Vega observations and the ATLAS9 Vega model 
were again used to define the zero points. 

The C96 comp osite spectrum represents both satellite 
dWoods et al.L [l996) and ground (Neckel & Labsl [19841) obser- 
vations below and above 4100 A. While its (b-y) is redder, m\ 
is considerably bluer, and c\ is in good agreement with our so- 
lar values. Not unexpectedly, the same conclusions about (b- 
y) and m\ hold for the N99 atlas, which is an update of the 
iNeckel & Labs ( 1984) measurements included in C96. However, 
Neckel (2003, hereafter N03) noticed a possible systematic er- 
ror in those absolute measurements and provides a simple ana- 



lytical formula for correction, after whic h (b-y) and m\ are in 
closer agreement with our solar values. Thuil lier et alJ (2004) 
published two composite solar reference spectra assembled us- 
ing space measurements during distinct solar activity levels; dif- 
ferences concern only the m\ and c\ indices in a negligible man- 
ner (a few millimagnitudes only), a nd therefore in Table 5 the 
averaged values of the two are given. IThuillier et al. (2004) pro- 
vide closer agreement than previous spectra in (b-y) and ni\, but 
their c\ is considerably bluer. 

Thus, while on average (b-y) and ;«i are in agreement with 
our solar values, there is the tendency for c\ to be considerably 
bluer. The /3 indices were computed following the same prescrip- 
tions as in the previous section and are systematically redder. 
Switching to the model flux of Vega solves most of the discrep- 
ancies for these two indices, with the synthetic values being now 
distributed at the red and blue sides of our solar colors (see the 
right panel of Fig. 5). It is not obvious why a model flux for Vega 
also brings absolute solar flux measurements into closer agree- 
ment with our empirical solar colors. While the measured spectra 
of Vega should be superior to any model of it, we can speculate 
about possible reasons for this not being true over the full wave- 
length range. The difference STIS vs. ATLAS9 amounts to few 
millimagnitudes in (b-y) and mi, but is considerable in c\ and 
/?, and those differences could reflect measurement problems in 
the u band and the issue already mentioned about the spectral 
resolution around the Balmer line. 

In the wavelength range of interest to us, the accuracy of the 
solar absolute fluxes is on the order of 5 - 3% from data ac- 
quired in space (C96; Thuilli er et all 12004) and probably lower 
for ground-based data measurements. The culprit could thus be 
inaccuracy in solar flux measurements rather than in Vega! For 
the sake of computing indices, we are not interested in the abso- 
lute flux scale, but rather in the relative accuracy of measuring 
the shape of the solar spectral energy distribution. This accuracy 
is rather difficult to assess, but the above differences provide an 
idea of the complexity associated with this kind of measurement, 
and the corresponding rather large uncertainties. In addition, the 
choice for Vega and the set of filters used also affect the analysis, 
introducing systematic errors. 

In this section, we have presented various issues concerning 
synthetic solar colors, which are often taken for granted when 
comparing observed vs. synthetic colors. While the comparison 
with a specific solar-flux measurement has limited significance, 
with appropriate choices for Vega, it is comforting that there are 
no large systematic trends between our solar colors and these ob- 
servations taken together. Thus, our solar colors can help consid- 
erably in establishing the necessary reliability of synthetic col- 
ors. 



6.3. Testing the effects of changes in T eS , iogg, and [Fe/H] 
on synthetic MARCS colors 

As discussed in the previous sections, our uvby-fi solar colors 
are useful for testing absolute solar spectra and the performance 
of model atmospheres, as well as checking the subtleties in the 
transformation from fluxes to synthetic colors. Since our solar 
twins and analogs span a range in T e s, log g, and [Fe/H], we 
can also use them to test how well model atmospheres predict 
the relative change of colors due to variations in atmospheric 
parameters. 

For this test, we used relative fluxes predicted by the 
MARCS models, normalizing the results to our "observed" so- 
lar colors. As shown in Figs. 1-4, where the relative variation 
in the MARCS (b-y), /3, mi, and c\ colors are shown by solid 



J. Melendez et al.: The colors of the Sun 



0.10 



0.05 



E 0.00 
< 



-0.05 - 



-0.10 



Colino et ol. (1996) ATLAS9 

Neckel (1999) MARCS 




Neckel (1999 + 2003) 
Thuillier et ol. (2004) 




(b-y) 



m. 



C1 



fi 



(b-y) 







Fig. 5. Comparison between solar and synthetic colors obtained from different spectra. Left panel: using Vega STIS observations 
to set the zero-points. Right panel: using ATLAS9 Vega model for the same purpose. Amag = our derived solar colors minus the 
synthetic ones. 



lines, there is a satisfactory overall relative agreement between 
MARCS models and both the solar twins and solar analogs. The 
only clear discrepancy is in the predicted metallicity variation 
in p (Fig. 2), which is much shallower than observed. We also 
checked whether the MARCS /3 colors computed by O09 solve 
this problem (Fig. 2), but although they do help to slightly al- 
leviate the discrepancy (the Onehag et al. /3 colors are slightly 
steeper than ours), they also do not match the low f3 indices ob- 
served below [Fe/H] ~ -0.2. It would be important to observe 
more metal-poor solar analogs to confirm whether we have de- 
tected a potential problem in the line formation of H/3. 

The strong variation in C[ with metallicity for G-type stars 
dTwarog et al., 2002, Fig. 4; see also O09) is reproduced well by 
the MARCS models (as already demonstrated by O09 for some- 
what cooler dwarfs), so it can potentially be used to estimate 
[Fe/H] in solar metallicity and metal-rich solar analogs. 

Our tests provide confidence in the relative variation in 
stellar fluxes with atmospheric parameters as predicted by the 
MARCS models, at least in the atmospheric parameter space 
studied here. 



7. The zero-point of the r eff scale 

Our solar (b-y) color can also be used to estimate how well dif- 
ferent T e g scales reproduce the solar value of T^ (5777 K). 
Different color- T^ relations in the literature have been em- 
ployed to obtain the T e ff for our (b - y) Q , and these values have 
then been subtracted from the solar effective temperature. These 
differences, shown in Table 6, are the zero-point corrections 
needed to place different T^ calibrations on t he same T e g scale 
as the Sun. As can be seen, the T e a scale by Casagrande et al.l 
(2010) seems the most accurate, with a negligible offset of only 
10(+10) K with respect to the Sun, which is much smaller than 

the offset (109 K) of our earlier T eS calibra tion (RM05b) . 

The widely used T e « (b-y) scale of Alo nso et al.l d 19961) 
is too cool by 137 K. Thus, the elemental chemical abun- 
dances determined using this sca l e (e.g . [israelian et al., 1998; 
Chen et all l2000t iReddv et all 120031: lAllende P rieto et all 



as much as 0.14 dex (e.g. for V from VI lines and for N from 
NH lines), albeit the impact on abundance ratios [X/Fe] should 
be smaller (and in a few cases mostly cancel) for lines of low to 
moderate excitation potential because of the compensating im- 
pact of the iron abundance obtained from Fel lines; however, for 
abundances obtained from lines of high excitation potential (e.g. 
CI, NI, Ol, PI, SI) the impact could be as high as 0.2 dex. Thus, 
inaccurate temperature scales may be one of the reasons why 
the small peculiarities in the solar chemical composition (M09; 
R09) have not been discovered in the past (Ramirez et al., 2010). 
Zero-point errors in T s $ also affect stellar ages inferred from 
isochrones. For exampl e, for an o ffset of 130 K and using Y 2 
isochrones (Demarqueet al., 2004), a systematic offset as high 
as 1.5-2 Gyr may result for a solar analog. 



8. The zero-point of the metallicity scale 



2004t iJonsell et all 120051) may be systematically incorrect by 



Arnadotti r et al.l (1201 Ol) made a critical evaluation of different 
metallicity calibrations for dwarf stars in the uvby system, find- 
ing zero-point offsets in many of these meta llicity scales, rang- 
ing from -0.17 dex dNordstrom et all 12004 to +0.33 (O09), 
although the latter is a purely theoretical calibration intended 
to test model atmospheres and not for application to real stars. 
After extensive testing for potential problems in different metal- 
licity calibrations, such as trends in the A[Fe/H] residuals 
(spectroscopic-photometric) with [Fe/H], (b - y),co, and mo, 
lArnadottir et al.l (1201 Ol) find that the calibrations by RM05a has 
the best overall performance, albeit with a small (0.04 dex) zero- 
point offset. 

The zero-point offset cannot be overstated, e specially when 
comparing the Sun to the stars (Gustafsson, 2008), because off- 
sets in the zero-point of the metallicity scale could make the Sun 
appear abnormal (see Haywood, 2008; Gustafsson et al., 2010; 
and references therein for a discussion about apparent anomalies 
in the solar metallicity). 

Our solar colors can be used to check the zero-point of dif- 
ferent metallicity calibrations, by computing [Fe/H] for our in- 
ferred solar colors and subtracting these metallicities from the 
solar metallicity ([Fe/H] = 0). Those differences, which are the 
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zero-point corrections needed to place different metallicity cali- 
brations on the same metallicity scale as the Sun, are shown in 
Table 7, along with the offsets found by Arnado ttir et alj {2010) 
for a broad range of colors and metallicities. 

Most metallicity calibrations do require a correction to place 
them on the same metallicity scale as the Sun. In particular, the 
RM05a metallicity calibration recommended by Arnadottir et al. 
(2010) to derive |Fe/H] from Stromgren photometry, needs a 
zero-point correction of +0.04 dex. This offset is identical to 
the offset found by Arnadotti r et alj d2010l) using a sample span- 
ning a much broader range in T e ff. We are currently revis- 
ing the metallici t y calibration employed in t he GCS survey 
dNordstrom et all l2004t iHolmberg et"aU 120071) and plan to as- 
sess the apparent anomalously high metallicity of the Sun with 
a new accurate metallicity scale (Casagrande et al., in prepara- 
tion). 



9. Conclusions 

New uvby-fi photometry has been presented for solar-twin can- 
didates observed at the SPM observatory. Comparisons with ex- 
isting Stromgren photometry shows that our data are in excellent 
agreement (to better than 0.01 mag) with previous observations. 

Using accurate spectroscopically derived stellar parameters, 
the uvby-[i photometry for the solar twins, and also for solar 
analogs covering a wider range in metallicities, the solar col- 
ors (b - y) = 0.4105+0.0015, m le = 0.2122+0.0018, c hG = 
0.3319+0.0054, and/? = 2.5915+0.0024 have been inferred. 

As discussed in the manuscript, our solar-twin data have pro- 
vided stringent constraints on absolute fluxes, the performance 
of model atmospheres, and the zero-points of temperature and 
metallicity scales. I n particular, we show that the widely used 
lAlonso et aD (l996) (b-y) calibr ation is too cool by ~ 140 K, 
while our new T e ff calibration (Casagrand e et all 1201 Oh has a 
negligible zero-point offset (10+10 K). Regarding model atmo- 
spheres, the Kurucz ATLAS9 solar model provides a closer fit of 
our solar colors than the MARCS 2008 solar model. The relative 
variation in colors with stellar parameters seem to be reproduced 
well by MARCS models, except for the metallicity variation in 
the fi index, which is flatter than observed in solar twins and so- 
lar analogs, thus suggesting that there are existing limitations on 
the modeling of Balmer lines. 

We are pursuing photometric observations of our solar-twin 
sample in other photometric systems, to perform similar evalu- 
ations of absolute fluxes, model atmospheres, and fundamental 
calibrations in astrophysics. 
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color 


A 


error 


B 


error 


C 


error 


D 


error 


a- (fit) 


solar colors 


b-y 

Cl 


0.4105 
0.2122 
0.3319 
2.5915 


0.0005 
0.0006 
0.0018 
0.0008 


-1.2737e-4 
-1.1405e-4 

4.5183e-5 


9.2777e-6 
1.1657e-5 

1.5840e-5 


0.05051 
-0.11223 
-0.03241 


0.00936 
0.02658 
0.01192 


0.049813 
0.125539 
0.124344 
0.034604 


0.004753 
0.005877 
0.016006 
0.007492 


0.0046 
0.0056 
0.0165 
0.0061 


0.4105+0.0015 
0.2122+0.0018 
0.3319+0.0054 
2.5915+0.0024 



Table 6. A T e fi needed to correct the zeropoint of the most recent 
r e ff scales 



AT tI 


reference 




(K) 






(b-y) 


+ 137 


Alonsoetal. (1996) 




+85 


Gratton et al. (1996) 




+118 


Blackwell & Lvnas-Grav 


(1998) 


+68 


Clem et al. (2004) 




+ 109 


RM05b 




+48 


Holmberg et al. (2007) 




+ 10 


Casagrande et al. (2010) 





p 



+ 130 Alonsoetal. (1996) 

Table 7. A [Fe/H] needed to correct the zero-point of different 
uvby metallicity scales 

A [Fe/H] A [Fe/H] reference 

This work Arnadottir et al. 2010 

01sen(1984) 

Schust er & Nissen (1989) 

Haywood (2002) 

Mart ell & Laughlin (2002) 

Martell & Smith (2004) 

RM05a 

Holmberg et al. (2007) 

O09 



+0.06,+0.05 


+0.11+0.34 


+0.05 


+0.06+0.16 


-0.04 


+0.00+0.18 


+0.07 


+0.05+0.13 


+0.07 


+0.06+0.21 


+0.04 


+0.04+0.14 


+0.09 


+0.08+0.16 


+0.37 


+0.33+0.30 
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Table 1. Photometry (V Johnson and uvby-jS) of 73 solar-twin 
candidates from San Pedro Martir Observatory 



HIP number 


HD/BD number 


V 


(b-y) 


mi 


C\ 


N v 


^ uvby 


P 


N /5 


348 


225194 


8.600 


0.407 


0.188 


0.307 


4 


4 


2.584 


4 


996 


804 


8.184 


0.403 


0.190 


0.340 


4 


4 


2.598 


4 


1411 


1327 


9.088 


0.441 


0.211 


0.342 


4 


4 


2.579 


4 


2894 


+48 0182 


8.653 


0.415 


0.208 


0.291 


4 


4 


2.587 


4 


4909 


6204 


8.508 


0.403 


0.206 


0.319 


4 


4 


2.597 


4 


5134 


6470 


8.965 


0.402 


0.187 


0.281 


4 


4 


2.581 


4 


6407 


8291 


8.612 


0.411 


0.205 


0.296 


4 


4 


2.584 


4 


7245 


9446 


8.377 


0.423 


0.218 


0.328 


4 


4 


2.596 


4 


8507 


11195 


8.890 


0.413 


0.207 


0.317 


4 


4 


2.595 


4 


8841 


-13 0347 


9.232 


0.423 


0.209 


0.301 


4 


4 


2.593 


4 


9349 


12264 


7.975 


0.407 


0.210 


0.324 


4 


4 


2.594 


4 


10710 


+ 18 0289 


8.924 


0.400 


0.191 


0.302 


5 


5 


2.594 


5 


11728 


15632 


8.041 


0.414 


0.219 


0.333 


5 


5 


2.590 


5 


11915 


16008 


8.608 


0.406 


0.211 


0.314 


4 


4 


2.594 


4 


14614 


19518 


7.858 


0.400 


0.198 


0.303 


5 


5 


2.591 


5 


18261 


24552 


7.976 


0.390 


0.195 


0.325 


5 


5 


2.604 


5 


25670 


36152 


8.285 


0.413 


0.219 


0.380 


6 


6 


2.599 


6 


28336 


40620 


8.982 


0.411 


0.209 


0.294 


6 


6 


2.578 


6 


35265 


56124 


6.936 


0.395 


0.207 


0.349 


7 


7 


2.594 


7 


36512 


59711 


7.732 


0.403 


0.210 


0.316 


5 


5 


2.588 


5 


38072 


63487 


9.206 


0.406 


0.215 


0.351 


6 


6 


2.601 


6 


41317 


71334 


7.808 


0.409 


0.218 


0.341 


5 


5 


2.588 


5 


41832 


71779 


8.125 


0.399 


0.199 


0.322 


7 


7 


2.580 


7 


42872 


74645 


9.276 


0.402 


0.198 


0.363 


4 


4 


2.594 


4 


44324 


77006 


7.942 


0.390 


0.190 


0.322 


6 


6 


2.587 


6 


44935 


78534 


8.722 


0.413 


0.212 


0.344 


5 


5 


2.590 


5 


44997 


78660 


8.346 


0.412 


0.221 


0.322 


5 


5 


2.588 


5 


47990 


84705 


8.704 


0.406 


0.218 


0.327 


5 


5 


2.591 


5 


49572 


+30 1962 


9.288 


0.406 


0.203 


0.346 


5 


5 


2.591 


5 


49756 


88072 


7.549 


0.403 


0.220 


0.326 


5 


5 


2.593 


5 


50826 


+ 17 2213 


9.129 


0.403 


0.189 


0.283 


5 


5 


2.581 


5 


51337 


90733 


8.905 


0.399 


0.203 


0.354 


5 


5 


2.595 


5 


52040 


91909 


9.194 


0.406 


0.207 


0.288 


5 


5 


2.573 


5 


52137 


92074 


8.638 


0.415 


0.218 


0.317 


5 


5 


2.584 


5 


55409 


98649 


8.010 


0.407 


0.217 


0.356 


5 


5 


2.585 


5 


55459 


98618 


7.658 


0.406 


0.206 


0.341 


5 


5 


2.607 


5 


58303 


103828 


8.443 


0.411 


0.210 


0.314 


5 


5 


2.642 


5 


59357 


105779 


8.665 


0.389 


0.189 


0.318 


5 


5 


2.581 


5 


60314 


107633 


8.783 


0.409 


0.213 


0.362 


4 


5 


2.592 


5 


60653 


108204 


8.737 


0.401 


0.181 


0.352 


5 


5 


2.579 


5 


63048 


112257 


7.808 


0.423 


0.220 


0.351 


5 


5 


2.582 


5 


64150 


114174 


6.795 


0.412 


0.234 


0.335 


5 


5 


2.584 


5 


64497 


114826 


8.943 


0.411 


0.215 


0.351 


5 


5 


2.598 


5 


64713 


115169 


9.268 


0.398 


0.211 


0.355 


5 


5 


2.579 


5 


64794 


115382 


8.433 


0.408 


0.201 


0.321 


3 


3 


2.584 


3 


64993 


115739 


8.899 


0.405 


0.213 


0.341 


3 


3 


2.588 


3 


65627 


+47 2060 


9.136 


0.408 


0.199 


0.308 


3 


3 


2.590 


3 


66885 


119205 


9.305 


0.403 


0.180 


0.304 


3 


3 


2.562 


3 


70394 


+29 2529 


9.568 


0.431 


0.226 


0.305 


4 


4 


2.580 


4 


73815 


133600 


8.198 


0.409 


0.212 


0.328 


4 


5 


2.600 


5 


74341 


134902 


8.861 


0.416 


0.214 


0.335 


3 


3 


2.588 


3 


74389 


134664 


7.781 


0.395 


0.214 


0.392 


3 


4 


2.590 


4 


75528 


+47 2225 


9.785 


0.427 


0.192 


0.262 


4 


4 


2.568 


4 


75923 


138159 


9.184 


0.414 


0.212 


0.306 


2 


2 


2.588 


2 


77883 


142331 


8.739 


0.422 


0.225 


0.345 


2 


2 


2.594 


2 


77936 


234267 


9.334 


0.412 


0.151 


0.292 


4 


4 


2.574 


4 


78028 


+37 2687 


8.640 


0.410 


0.185 


0.307 


4 


4 


2.590 


4 


78680 


144270 


8.198 


0.404 


0.182 


0.299 


4 


4 


2.594 


4 


79186 


145514 


8.316 


0.427 


0.185 


0.322 


4 


4 


2.583 


4 


79304 


145478 


8.684 


0.412 


0.192 


0.359 


4 


4 


2.594 


4 


79672 


146233 


5.494 


0.400 


0.219 


0.356 


4 


4 


2.603 


4 


81512 


+45 2434 


9.228 


0.425 


0.189 


0.323 


4 


4 


2.581 


4 


85285 


157691 


8.388 


0.406 


0.164 


0.312 


4 


4 


2.583 


4 


88194 


164595 


7.070 


0.415 


0.195 


0.336 


4 


4 


2.591 


4 


88427 


+35 3136 


9.331 


0.418 


0.179 


0.314 


4 


4 


2.590 


4 


100963 


195034 


7.091 


0.408 


0.199 


0.327 


4 


4 


2.586 


4 


102152 


197027 


9.193 


0.417 


0.210 


0.355 


4 


4 


2.595 


4 


103025 


+ 14 4456 


8.719 


0.414 


0.190 


0.353 


4 


4 


2.587 


5 


104504 


201422 


8.550 


0.396 


0.191 


0.305 


4 


4 


2.589 


5 


108708 


209096 


8.943 


0.415 


0.217 


0.353 


4 


4 


2.595 


4 


108996 


209562 


8.894 


0.406 


0.209 


0.345 


4 


4 


2.598 


4 


109931 


+24 4563 


8.941 


0.423 


0.205 


0.350 


4 


4 


2.588 


4 


118159 


224448 


9.003 


0.400 


0.195 


0.295 


3 


4 


2.589 


4 
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Table 2. Literature photometry (V J ohnson and uvby-fi) of 
solar-twin candidates and solar analogs (M ermiHiod et al.Lll997t 
Hauck & Mermilliod, 1998b 



HIP number HD number 



V (b-y) 



ni] 



C] 



N v 



P N« 



996 


804 


8.191 


0.395 


0.196 


0.355 


2 






11728 


15632 


8.035 


0.414 


0.221 


0.335 


3 






22263 


30495 


5.489 


0.398 


0.213 


0.321 


15 


2.601 


11 


29525 


42807 


6.440 


0.415 


0.228 


0.292 


13 


2.593 


10 


30502 


45346 


8.660 


0.411 


0.202 


0.331 


4 






36512 


59711 


7.742 


0.406 


0.204 


0.308 


2 






38228 


63433 


6.891 


0.423 


0.232 


0.313 


1 






41317 


71334 


7.812 


0.412 


0.210 


0.327 


4 






42438 


72905 


5.651 


0.396 


0.197 


0.293 


4 


2.596 


4 


43190 


75288 


8.515 


0.423 


0.234 


0.345 


3 






44713 


78429 


7.303 


0.414 


0.227 


0.348 


3 


2.600 


19 


44997 


78660 


8.345 


0.410 


0.216 


0.321 


2 






49756 


88072 


7.544 


0.407 


0.211 


0.335 


4 






55409 


98649 


8.004 


0.405 


0.227 


0.323 


1 


2.588 


1 


55459 


98618 


7.658 


0.411 


0.198 


0.347 


2 






56948 


101364 


8.673 


0.410 


0.212 


0.313 


1 






64150 


114174 


6.791 


0.419 


0.231 


0.334 


15 


2.592 


9 


77052 


140538 


5.869 


0.424 


0.231 


0.334 


5 






79672 


146233 


5.496 


0.401 


0.217 


0.343 


11 


2.595 


3 


85042 


157347 


6.287 


0.422 


0.225 


0.346 


4 






100963 


195034 


7.090 


0.408 


0.204 


0.319 


1 






102152 


197027 


9.179 


0.421 


0.202 


0.323 


1 






109110 


209779 


7.581 


0.415 


0.216 


0.318 


2 






1499 


1461 


6.468 


0.421 


0.238 


0.363 


15 


2.597 


13 


15457 


20630 


4.850 


0.419 


0.235 


0.307 


54 


2.595 


10 


53721 


95128 


5.037 


0.391 


0.202 


0.343 


37 


2.606 


13 


59610 


106252 


7.425 


0.390 


0.187 


0.341 


3 






60081 


107148 


8.021 


0.424 


0.244 


0.394 


2 






62175 


110869 


8.008 


0.412 


0.215 


0.358 


3 






79578 


145825 


6.563 


0.395 


0.228 


0.328 


1 


2.597 


1 


80337 


147513 


5.373 


0.399 


0.199 


0.323 


5 


2.609 


1 


96402 


184768 


7.556 


0.426 


0.214 


0.343 


6 


2.590 


6 


96895 


186408 


5.979 


0.410 


0.212 


0.370 


64 


2.606 


39 


96901 


186427 


6.234 


0.416 


0.223 


0.352 


63 


2.597 


39 


113357 


217014 


5.456 


0.416 


0.232 


0.371 


143 


2.603 


13 
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HIP number 7^ log g [Fe/H] reference 



348 


5777 


4.41 


-0.13 


R09 


996 


5860 


4.38 


0.00 


R09 


1499 


5756 


4.37 


0.19 


R09+VF05+LH06+T07+S08 


2894 


5820 


4.54 


-0.03 


R09 


4909 


5836 


4.44 


0.02 


R09 


5134 


5779 


4.49 


-0.19 


R09 


6407 


5787 


4.47 


-0.09 


R09 


8507 


5720 


4.44 


-0.08 


R09 


8841 


5676 


4.50 


-0.12 


R09 


9349 


5825 


4.49 


0.01 


R09 


10710 


5817 


4.39 


-0.13 


R09 


11728 


5738 


4.37 


0.05 


R09+T07 


11915 


5793 


4.45 


-0.05 


R09 


14614 


5803 


4.47 


-0.10 


R09+T07+B10 


15457 


5771 


4.56 


0.08 


B10+VF05 


18261 


5891 


4.44 


0.00 


R09+T07 


22263 


5826 


4.54 


0.00 


B10+VF05 


25670 


5755 


4.38 


0.07 


R09+T07 


28336 


5713 


4.53 


-0.17 


R09 


29525 


5715 


4.41 


0.00 


BIO 


30502 


5745 


4.47 


-0.01 


M09 


36512 


5740 


4.50 


-0.09 


M09+T07+S08 


38072 


5839 


4.53 


0.06 


R09 


38228 


5693 


4.52 


0.01 


R09+VF05 


41317 


5724 


4.46 


-0.04 


M09+VF05+S08 


42438 


5864 


4.46 


-0.05 


R09 


43190 


5775 


4.37 


0.12 


M09 


44324 


5934 


4.51 


-0.02 


R09+T07 


44713 


5784 


4.36 


0.10 


B10+VF05+S08 


44935 


5800 


4.41 


0.07 


M09 


44997 


5773 


4.53 


0.03 


M09+T07 


49572 


5831 


4.33 


0.01 


R09 


49756 


5804 


4.45 


0.04 


R09+VF05+T07 


50826 


5725 


4.47 


-0.28 


M09 


52040 


5785 


4.51 


-0.09 


R09 


52137 


5842 


4.56 


0.07 


R09 


53721 


5916 


4.48 


0.03 


B10+VF05 


55409 


5760 


4.52 


-0.01 


M09 


55459 


5838 


4.42 


0.04 


R09+VF05+M06+MR07+T07 


56948 


5795 


4.45 


0.02 


R09+MR07+T09+M10 


59357 


5810 


4.45 


-0.24 


M09 


59610 


5899 


4.34 


-0.03 


B10+VF05 


60081 


5811 


4.38 


0.32 


M09+VF05+S08 


60314 


5874 


4.52 


0.11 


R09 


60653 


5725 


4.38 


-0.29 


M09 


62175 


5849 


4.43 


0.14 


R09+T07 


64150 


5755 


4.39 


0.06 


R09+VF05+T07 


64713 


5815 


4.52 


-0.01 


M09 


64794 


5743 


4.33 


-0.10 


R09 


64993 


5875 


4.56 


0.09 


M09 


66885 


5685 


4.48 


-0.38 


M09 


73815 


5803 


4.34 


0.02 


MR07+R09 


74341 


5853 


4.51 


0.09 


R09 


74389 


5859 


4.48 


0.11 


M09+S08 


75923 


5775 


4.56 


-0.02 


M09 


77052 


5697 


4.54 


0.04 


B10+VF05 


77883 


5695 


4.39 


0.04 


M09 


79578 


5860 


4.53 


0.07 


M09+VF05 


79672 


5822 


4.45 


0.05 


M09+VF05+M06+MR07+T07+T09+M09 


80337 


5881 


4.53 


0.03 


B10+VF05+S08 


85042 


5692 


4.39 


0.04 


B10+VF05+S08 


85285 


5730 


4.43 


-0.39 


M09 


88194 


5735 


4.40 


-0.07 


R09+VF05+T07 


88427 


5810 


4.42 


-0.16 


R09 


96402 


5713 


4.33 


-0.03 


B10+T07 


96895 


5808 


4.33 


0.10 


R09+VF05+LH06 


96901 


5737 


4.34 


0.06 


R09+ VF05 +LH06+T07 


100963 


5802 


4.45 


0.01 


R09+T07+T09 


102152 


5737 


4.35 


-0.01 


R09+M10b 


104504 


5836 


4.50 


-0.16 


R09 


108708 


5875 


4.51 


0.15 


R09 


108996 


5838 


4.50 


0.06 


R09 


109110 


5817 


4.46 


0.06 


B10+VF05+T07 


109931 


5739 


4.29 


0.04 


R09 


113357 


5803 


4.38 


0.22 


R09+VF05+LH06 
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Table 5. Solar colors and predictions from models atmospheres. 
Our synthetic colors (This work) were computed using the 
Bessell (2005) filters, with zero points based on both the STIS 
observed and ATLAS9 synthetic spectrum of Vega 



(H) 


m\ 


C\ 


fi 


reference 


Solar colors (this work) 


0.4105 (±0.0015) 


0.2122 (±0.0018) 


0.3319 (±0.0054) 


2.5915 (±0.0024) 




Solar spectra 


0.406,0.390(±0.004), (b-y) from/3 






2.5955(±0.0024) 


Olsen (1976), using asteroids 


0.425 


0.172 


0.323 


2.625 


This work, C96 [STIS] 


0.421 


0.178 


0.350 


2.595 


This work, C96 [ATLAS9] 


0.424 


0.188 


0.306 


2.632 


This work, N99 [STIS] 


0.420 


0.194 


0.333 


2.602 


This work, N99 [ATLAS9] 


0.400 


0.195 


0.305 


2.631 


This work, N99+N03 [STIS] 


0.396 


0.201 


0.332 


2.601 


This work, N99+N03 [ATLAS9] 


0.402 


0.225 


0.287 


2.625 


This work, Thuillier (2004) [STIS] 


0.398 


0.231 


0.314 


2.595 


This work, Thuillier (2004) [ATLAS9] 


Kurucz models 


0.371 


0.214 


0.243 




Relvea& Kurucz (1978) 


0.414 








Kurucz (1991) 


0.355, 0.388 


0.235 


0.359 


2.618, 2.660 

2.581 


Lester et al. (1986) 

Smallev & Dworetskv (1995) 


0.393(CM), 0.400(noOV), 0.414(OV) 


0.278 (CM) 


0.339 (CM) 




Smallev & Kuoka (1997) 


0.397(noOVER), 0.410 (OVER) 






2.590 


Castelli et al. (1997) 
Castelli & Kurucz (2006) 


0.410 
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